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Abstract. The paradigm for singular continuous spectra in symbolic dynamics and in 
mathematical diffraction is provided by the Thue-Morse chain, in its realisation as a bi- 
nary sequence with values in {±1}. We revisit this example and derive a functional equation 
together with an explicit form of the corresponding singular continuous diffraction measure, 
which is related to the known representation as a Riesz product. 

The Thue-Morse chain can be defined via the primitive substitution rule 

1 i — y 11 
g: I-Il 

on the two-letter alphabet {1, 1}; see [HQI)] for background. Let v = v v 1 v 2 ... be the unique 
one-sided fixed point of g with v Q = 1. This infinite word satisfies 

(!) v 2n = v n and v 2n+l = 

for all n £ N , where 1 = 1. Consider the bi- infinite word w, defined as 

{v n , for n > 0, 

v- n -i, for n < 0, 

which is the unique reflection-symmetric fixed point of g 2 with (admissible) seed w_ 1 w = 11. 
It is cube-free and thus aperiodic. If S denotes the (two-sided) shift, defined by (Sw) n = 
w n+ i, the Thue-Morse hull is the compact space X = {S m w \ m £ Z}, where the closure 
is taken in the obvious product topology. Note that X coincides with the (discrete) local 
indistinguishability (or LI) class of w, and that the topological dynamical system (X, S) is 
strictly ergodic; see [10] for details. 

Here, we consider the weighted Dirac comb 

(2) u = w n S n , 

where 5 X denotes the normalised Dirac measure on M, located at x, and 1 = —1. In particular, 
u is a translation bounded, signed measure on M. The corresponding autocorrelation (or 
Patterson) measure 7, obtained by a volume-averaged convolution [61 [2] of u with its reflected 
counterpart, reads 

(3) 7 = Vh Yl ^( m )^' 
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with the coefficients rj(m) obtained as the limits 



1 1 N 
ri(m) = lim — — > w h w f = lim — — > w n w r , 



N^oo 2N + 1 ^ AT-^oo 2iV + _ 

-N<k,i<N n=-N 
k—£=m 

Note that these limits exist due to unique ergodicity, compare |10| . with r/(0) = 1. The 
autocorrelation coefficients satisfy rj(—m) = 77 (m) and, due to ([1]), the recursions 

(4) rj (2m) = rj (m) and rj(2m + 1) = — — (r](yri) + 77(771 + 1)) , 

valid for all m > 0. In particular, 77(1) = — |, and all coefficients are uniquely specified. 

By construction, r/ is a positive definite function on Z, wherefore the Herglotz-Bochner 
theorem [SI Thm. 1.7.6] guarantees the existence of a finite positive measure v on [0, 1) with 

(5) n(m) = f\ 2 ™ my du(y). 



The diffraction measure 7 of the weighted Thue-Morse comb u of ([2|) is the Fourier transform 
of 7. An elementary calculation shows that 



7 = v * 5% and v = 7 



[(VO- 



ID, this formulation, the spectral properties of 7 follow immediately from those of v. In 
particular, v has a unique decomposition 

v = U PP + ^sc + ^ac 

into its pure point, singular continuous and absolutely continuous parts (relative to Lebesgue 
measure A); see |12} Thms. 1.13 and 1.14]. 

It was first shown by Kakutani |7J that v is a purely singular continuous measure, so 
v = v sc . Let us briefly adapt this to our setting. By Wiener's criterion P Cor. 7.11], v pp = 
is equivalent to limAr^oo 2jV 1 +1 E(N) = 0, where 

N 



E(N) := H m )) ■ 



m=-N 

When N > 1, we can use the recursion relations @ to derive U(4N) < § 27(2JV). With 
a = log 2 (3/2) < 1, one then obtains the estimate j^U(N) = 0{\/N 1 ~ a ). At this point, we 
know that u = u sc + f ac , and define the distribution function 

(6) F(x):=u([0,x}), 

which is a continuous function of bounded variation on [0,1]. It satisfies F(0) = and 
F(x) + F(l — x) = 1, the latter as a consequence of the symmetry of v. 

Following [7] and viewing v as a Lebesgue-Stieltjes measure [HI Ch. X] with distribution 
function F, the two recursion relations (j3|) imply the identities 

dF(|) + dF(2±i) = dF(x), 

<LF(§) - dF(2±i) = -cos(7ra;)dF(x), 
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for all x G [0,1]. The left-hand sides are obtained from the corresponding sides of the 
recursions by a change of variables, followed by a split of the new integration region into 
two parts. The actual equality of the measures follows because we obtain equality of the 
integrals over arbitrary trigonometric polynomials, whence the Fourier uniqueness theorem 
[8j Thm. 1.2.7] applies. The relations ([7]) must also hold for the two continuous components 
of v separately, because v sc _!_ v ac by [9, Thm. VII. 2. 4]. 

Writing u ac = gX with g G -^ 1 ([0, 1]), the identities ([7]) result in 

l{9(l)+9(^ k ))=9(x) and \ ( 5 (f) - g{^)) = - cos(Trx) g(x), 

this time for almost all x G [0,1]. Defining now r] ac (m) = J Q e 2mmx g(x) dx, we have 
Vac(~ m ) = Vac( m ) an d inherit a set of recursions identical to however with rj ac (0) = as 
a result of the Riemann-Lebesgue lemma [8l Thm. 1.2.8], and hence r] ac (m) = for all meZ. 
By the Fourier uniqueness theorem [8j Thm. 1.2.7], this implies g{x) = almost everywhere, 
and hence v ac = 0. Since v itself cannot be the zero measure, we have v = v sc ^ 0, and 
the Thue-Morse diffraction measure 7 is a purely singular continuous, Z-periodic, positive 
measure. It is the same for all Dirac combs of members of the hull X. 

The remainder of this note is concerned with an explicit calculation of u, via its distribution 
function F, which does not seem to have been calculated in the literature so far. Adding the 
two equations for dF from ([7J), followed by integration, results in the functional equation 

1 f 2x 

(8) F(x) = - J (l-cos(7ry))dF(y) 

for the continuous function F, valid for x 6 [0, 5]. Since F is a continuous non-decreasing 
function on [0, 1], the difference f(x) = F(x) — x defines a continuous function / of bounded 
variation, with /(0) = and f(x) + f(l — x) = for all x £ [0, 1]. As such, / possesses a 
uniformly converging Fourier series of the form 

f( x ) = b m sm(2irmx) , 

m>l 

with b m = 2 Jq 1 sin(27rmx) f(x) dx. From Eqs. ([5]) and ([!]), one finds rj(m) = irmb m , so that 

(9) Fix) = x + V sin(27rmx) 

^-^ ran 

m>l 

is a uniformly converging Fourier series representation of the Thue-Morse distribution func- 
tion. This function satisfies the symmetry relation F(x) + -F(l — x) = 1 for all x G [0, 1], and 
is a solution of the integral equation ([8]). 

To interpret ([8]) in a wider setting, let us introduce the space 

(10) d = |g g c([o, i],r) I G ^ + °^-°:^ 

of continuous, non-decreasing distribution functions on [0, 1] with the required symmetry, and 
define the mapping ^ : D — > D via G 1— ► ^(G) with 



(*(G))C> 



'i/ 2x (l-cos(vry))dG(y), if0<x<i, 
1- (&(G))(l-x), ifi<x<l. 
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One can show that F of Q is the only solution of $(F) = F within D, because <P is a weak 
contraction with respect to the metric d(G,H) := V(G — H), where V denotes the total 
variation of continuous functions on [0, 1]; see [HI Ch. X] for background. 

For a numerical high-precision computation of F, we employ a Volterra-type iteration 
within D, which is superior to using © directly. Starting from Fq(x) = x, where Fq G D, we 
set F n+ i = ^{F n ) for re > 0. This defines the (uniformly converging) sequence of distribution 
functions in D given by 

2 n -l (n) 

(11) F n (x) = x + } — sin(27rmx). 

mir 

m=l 

Here, the coefficients Sm satisfy the recursions 

(n+l) _ („) , (n+1) _ _\(Sn) , M \ 

c 2m — L m dlm c 2m+l ~ 2^ m rn+1! 

for n > 1 and < m < 2 n - 1, together with = 0, c[ 1} = -\ and c^ +1) = -\{\ + cj n) ). 
They approach the autocorrelation coefficients 77(771) via lirrin^oo cffl = 7/(m) for arbitrary 

777 E N. 

By construction, the distribution functions F n represent absolutely continuous measures 
(although their limit does not). Writing dF n (x) = f n (x)dx, one finds f (x) = 1, and the 
functional equation (|8|) induces the recursion 

f n+1 (x) = (l-cos(2vrx))/ n (2x) = 2(sin(vrx)) 2 /n(2x) 

for n > 0. This gives the well-known explicit representation 

n— 1 n— 1 

(12) / n (x) = J] C 1 - cos(2 <; + 1 ^)) = T H (sin(2^x)) 2 

£=0 1=0 

of the Thue-Morse measure as a Riesz product; compare [HI Sec. 1.4.2]. 

The Volterra iteration leads to a sequence (F n ) ng N of continuous distribution functions 
that converge uniformly (on [0, 1]) to F, which is continuous as well. The latter is shown in 
Figure [T]and resembles the classic middle-thirds Cantor measure in various aspect, though it 
has full support (see below). Note that the corresponding sequence (d-F n ) ng N o of absolutely 
continuous measures is only vaguely convergent, with the limit being purely singular continu- 
ous. Therefore, it is somewhat misleading to show a density for the Thue-Morse measure, as 
is often found in the literature. Still, it may be instructive to inspect the sequence of densities 
f n in order to get some intuition on the singular nature of the Thue-Morse measure, or to 
study some of its scaling properties; see [5] and references therein. 

On initial inspection, the distribution function F seems to have a plateau around x = |, 
with exact value ^. More generally, as suggested by the Riesz products (fl~2|) . one might expect 
a plateau around any x G ({0, 1} U {m/2 k j k € N, m G N odd}) PI [0, 1], because the densities 
f n vanish at x for all sufficiently large 77, with the order of this zero linearly increasing with 77. 
However, these potential gaps are all closed (see below). The corresponding values of F can 
be calculated with the series expansion Q. Unlike the situation of the gap labelling theorem 
for one-dimensional Schrodinger operators [3], where one knows the values of the integrated 
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Figure 1 . Distribution function F of the purely singular continuous diffrac- 
tion measure of the Thue- Morse chain in the interval [0, 1]. Note that F is a 
strictly increasing function, despite the appearance of plateau-like regions. 

density of states (IDOS) on the (always non-overlapping) plateaux but not their positions, 
we know the possible locations of the plateaux, but cannot see a topological constraint for 
the corresponding values of the distribution function (which can be determined from ([9])). 

It is interesting to note that the set of potential plateau locations coincides with the set of 
potential (but in our case extinct) Bragg peak positions, so the (extinct) Bragg peaks appear 
to 'repel' the continuous diffraction spectrum. Nevertheless, one has 

supp( dF) = supp(^) = [0, 1] 

by [U Prop. 28], which also implies that F is a strictly increasing function. In particular, 
as one can immediately see from ([8]), F(x) = forces F(2x) = 0, which (when repeated) 
contradicts F(l/2) = 1/2 unless x = 0. This shows that there is no gap around (and none 
around 1 by symmetry). The general argument uses that v is a regular Borel measure; see 
[1] and references therein for details. 
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The methods used above can also be applied to other substitutions of constant length that 
fail to have pure point spectrum. This can be decided on the basis of Dekking's criterion; 
see Sec. 6] for details. Although one still has to check Wiener's criterion and to find the 
analogue of the functional equation ([7]), this approach seems worth pursuing. 
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